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Locally rotationally symmetric (LRS) Bianchi Type I cosmological models are examined 

in the presence of dynamically anisotropic dark energy and perfect fluid. We assume that 

^ , 1 the dark energy (DE) is minimally interacting, has dynamical energy density, anisotropic 

equation of state parameter (EoS). The conservation of the energy-momentum tensor of 

$H , the DE is assumed to consist of two separately additive conserved parts. A special law is 

OXJ. assumed for the deviation from isotropic EoS, which is consistent with the assumption on 

the conservation of the energy-momentum tensor of the DE. Exact solutions of Einstein's 

^s^ , field equations are obtained by assuming a special law of variation for the mean Hubble 

^ ' parameter, which yields a constant value of the deceleration parameter. Geometrical and 

f^^ ' kinematic properties of the models and the behaviour of the anisotropy of the dark energy 

VO ■ has been carried out. The models give dynamically anisotropic expansion history for the 

OO ' universe that allows to fine tune the isotropization of the Bianchi metric, hence the CMB 

\l I anisotropy. 

00 ! 1 Introduction 

O 

One of the most successful attempt to resolve the problems of standard Big Bang cosmology 
• r^ . such as homogeneity, isotropy, and flatness of the universe is the inflationary paradigm, charac- 

^ ' terised by an epoch of accelerated expansion, " inflation" , in the very early universe [TJ [21 [31 H] . 

^ . During the inflationary epoch, quantum fluctuations are highly amplified, their wavelengths are 

stretched to outside the Hubble horizon and inevitably, superhorizon fluctuations are gener- 
ated. These fluctuations become classical after crossing the event horizon and are coherent on 
what appear to be superhorizon scales at decoupling [H [SI [71 [51 [HI [IHl [HI [HI [131 E] • After 
the end of inflation, they re-enter the horizon, and seed the matter and the radiation fluc- 
tuations observed in the universe. These primordial fluctuations are Gaussian, adiabatic and 
nearly scale-invariant in the majority of inflation models and even a slight deviation from these 
properties can strongly constrain the assumptions in cosmological models [13l [14] . Since the 
linearity of the cosmic microwave background (CMB) anisotropy preserves the basic properties 
of these primordial fluctuations, the CMB radiation anisotropy is a promising tool for testing 
these properties [13]. However, it is known that the observed quadrupole (/ — 2, which is the 
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largest scale that can be observed) amplitude has a lower value than the quadrupole expected 
from a best-fit A-dominated cold dark matter (ACDM) standard model to the entire power 
spectrum since the first data of the differential microwave radiometer (COBE/DMR) appeared 
in 1992 [121 [TB]. This anomaly was confirmed with the high resolution data provided by the 
first year (2003) and 3-year (2006) Wilkinson Microwave Anisotropy Probe (WMAP, WMAP3) 
[TTIIIH]- Most recently, the 5-year WMAP data did not improve on the quadrupole, it still seems 
to be outside of the fit [19l [20] . Today, this low value of the quadrupole seems inescapable and 
the main issue now is to elucidate possible explanations [H]. The proposed four possible candi- 
dates are that this anomaly is due to a systematic error, a pure statistical fiuke, astrophysical 
(i.e. unexpected foreground) or cosmological reasons (e.g. nontrivial spatial geometry of the 
universe) [3TJ[52]. Although the correct explanation is still unknown, the first three candidates 
seem unlikely [20l [211 [22l [23l [24l [25] . If it is indeed due to cosmological reasons, since the low- 
est multipoles represent especially the scale of the horizon at approximately dark energy (DE) 
domination, it seems natural to associate the low value of the quadrupole with the nature of 
the DE. 

It is known that Bianchi universe anisotropics give rise to CMB anisotropics depending 
on the model type [56]. Recently, Campanelli et al. showed that allowing large-scale spatial 
geometry of the universe to be plane-symmetric with eccentricity (regardless of the origin) at 
decoupling of order 10~^ can bring the quadrupole amplitude in accordance with observations 
without affecting higher multipoles of the angular power spectrum of the temperature anisotropy 
[271 [55] ■ Although they are controversial, it should be mentioned here that there are also some 
independent indications of a symmetry axis in the large-scale geometry of the universe, coming 
from the analysis of spiral galaxies in the Sloan Digital Sky Survey (SDSS) [29] and from the 
analysis of polarization of electromagnetic radiation propagating over cosmological distances 
[30] [31] [32] [33l [34l 135] [36] [37] . A spatially ellipsoidal geometry of the universe can be described 
with Bianchi type metrics. However, Bianchi type I, V, VII models isotropize at late times even 
for ordinary matter, and the possible anisotropy of the Bianchi metrics necessarily die away 
during the inflationary era [26l [38]. In fact this isotropization of the Bianchi metrics is due 
to the implicit assumption that the DE is isotropic in nature. If the implicit assumption that 
the pressure of the DE is direction independent is relaxed, the isotropization of the Bianchi 
metrics can be fine tuned to generate arbitrary ellipsoidality (eccentricity). Therefore, the CMB 
anisotropy can also be fine tuned, since the Bianchi universe anisotropics determine the CMB 
anisotropics. The price of this property of DE is a violation of the null energy condition (NEC) 
since the DE crosses the Phantom Divide Line (PDL) , in particular depending on the direction 
[39]. 

However, it has been known since the 1980s that such energy components might occur and 
their role as possible DE candidates was raised by Caldwell at the end 1990s (see [40] for further 
references). In theory, despite the observational constraints, extensions of general relativity 
are the prime candidate class of theories consistent with PDL crossing [41]. On the other 
hand, while the current cosmological data from SNIa (Supernova Legacy Survey, Gold sample 
of Hubble Space Telescope) ^\^, CMB (WMAP, BOOMERANG) [H [1^ and large scale 
structure (SDSS) [JS] data rule out the w ^ —1, they mildly favor dynamically evolving DE 
crossing the PDL (see jJOJ [S] ; O [IH] for theoretical and observational status of crossing the 
PDL). 

Recently, Rodrigues [39] and Koivisto & Mota |49l [50] have investigated cosmological models 
with anisotropic equation of state (EoS). Rodrigues has constructed a Bianchi type-I ACDM 



cosmological model with a DE component which is non-dynamical but wields anisotropic vacuum 
pressure in two ways: i) by implementing of anisotropic vacuum pressure consistent with energy- 
momentum tensor conservation; ii) by implementing a Poisson structure deformation between 
canonical momenta such that rescaling of the scale factor is not violated [39] • He suggests to fine 
tune the DE so as to not wipe out the anisotropic imprints in the inflationary epoch. On the other 
hand, Koivisto & Mota have proposed a different approach to resolve CMB anisotropy problem; 
even if the CMB formed isotropically at early time, it could be distorted by the direction 
dependent acceleration of the later universe in such a way that it appears to us anomalous at 
the largest scales. They have investigated a cosmological model containing a DE component 
which has a non-dynamical anisotropic EoS and interacts with the perfect fluid component. They 
have also suggested that cosmological models with anisotropic EoS can explain the quadrupole 
problem and can be tested by SNIa data [49l [50] . 

In reference [51] Mota et al. have concluded that even though a perfect fluid representation 
might ultimately turn out to be a phenomenologically sufficient description of all the observa- 
tional consequences of DE, imperfectness in DE cannot be excluded J51]. Although there is 
compelling evidence that the expansion of universe is speeding up, we are far from understand- 
ing of the nature of the DE which is thought be the reason for this behaviour [JUIIS^FSH] . Hence, 
we should examine models with anisotropic dark energy, in order to determine what possible 
new physical consequences they might give rise to, and if for no other reason than to rule such 
models out. 

The above discussions lead us to examine the physical behaviours of the two-fluid, par- 
ticularly ACDM, cosmological models with an anisotropic DE component. Models of DE are 
conveniently characterized by the EoS parameter w — p/p which is not necessarily constant, 
where p is the energy density and p is the pressure [40]. However, while energy density is a 
scalar quantity, pressure is a vectorial quantity, and consequently the EoS parameter of DE 
may be determined separately on each spatial axis in a consistent way with the conservation of 
energy-momentum tensor. Hence, we consider a phenomenological parametrization of minimally 
interacting DE in terms of its time-dependent deviation-free equation of state parameter w^'^°' {t) 
and deviation parameters [5{t),j{t),j{t)]. Since such a parametrization yields an anisotropic 
expansion which is not compatible with the Robertson- Walker (RW) metric, in section 2 we have 
used the locally rotationally symmetric Bianchi Type I (LRS Bianchi-I) metric which generalizes 
the flat RW metric in an axially symmetric way and is compatible with our parametrization. The 
other component has assumed to be a perfect fluid (dark matter or ordinary matter). With our 
assumptions, the conservation of the energy-momentum tensor implies dynamically anisotropic 
DE except for very special solutions. We have obtained exact solutions for the equations by as- 
suming a special dynamic for the anisotropy of the dark the energy and a special law of variation 
for the mean Hubble parameter in Bianchi metrics, which yields a constant value of deceleration 
parameter [SH [551 \EM ■ The assumption on the mean Hubble parameter allows us to determine 
the scale factors exactly, as well as to examine the behaviour of the anisotropy of DE and other 
cosmological parameters of such a universe. 



2 Model and field equations 

The spatially homogenous, anisotropic and LRS Bianchi-I space-time is described by the line 
element 

ds^ = dt^ - A{t)^dx^ - B{tf{dy^ + dz^) (1) 

where A{t) and B{t) are the scale factors (metric tensors) and functions of the cosmic time t. 
In natural units (SttG = 1 and c = 1), the field equations, in the case of a mixture of the perfect 
fluid and the anisotropic DE components, are 

n — R /?„ — _7^(m) _ 7^(do) /2\ 

with 

T^-")/ = diagip^-"), -pf-"), -p^-"), -p^"")] (3) 

= diag[l, -w^'^\-w^'^\ -w;('")]p('") 

and 

r(^^)/ = diag[p('i°), -pj^^\-py(''^\-pj''^^] (4) 

- diag[l, -wj^^'l-wy^^^l -«;,(^^)]p(^^) 
= diag[l, -(zi;(^^) + (5), -(M;('i°) + 7), -(«;(^^) + 7)]p(^^) 

where g^i.u'^u'^ ~ 1; u'^ = (1,0,0,0) is the four-velocity vector; R^^ is the Ricci tensor; R is 
the Ricci scalar; p^™'' and p'-'^°^ are the energy densities of the perfect fluid and DE components 
respectively; w'™) is the EoS parameter of the perfect fluid and we will call u'^'^'^' the deviation- 
free EoS parameter of the DE. Here, the perfect fluid represents the ordinary matter or cold dark 
matter, thus w^™) > 0. Wx^'^'^\ Wy^'^'^'^ and Wz'^'^°^ are the directional EoS parameters of the DE 
on X, y and z axes respectively. Wy'^'*°-' and Wz'''^°^ are set to be equal which is convenient with 
the metric given in (1). 5 and 7 are deviations from the deviation- free EoS parameter (hence 
the deviation-free pressure) of the DE respectively on x axis and y and z axes. In a comoving 
coordinate system, Einstein's field equations (2), for the anisotropic LRS Bianchi-I space-time 
(1), in case of (3) and (4), read as 

— + 2-— -»('") + fl('i°) (^\ 

|J + 2| = -«;(-)p(-)-(^«('i°)+<5)p('^°^ (6) 

§ + §j + J = --^'"V^-' - (-('^"^ +7)p('^°\ (7) 

where the over dot denotes derivation with respect to the cosmic time t. We have the following 
equation from the Bianchi identity, 

G^^,=t(-)% + t(^^)^':. = 0, (8) 



which yields 



p(-) + (l + ,^(-)jp(-) - + 2- (9) 



+p(^^) + (l + u;('i°)) p^"^^ ( ^ + 2f ) + ''*'°^ (4 + ^^ 






This equation, which is hnearly dependent to Einstein field equations, also represents the con- 
servation of the total energy momentum tensor. 

3 Solution of the field equations 

We have initially 8 variables [A, B, p^™-', w^™-*, p^'^'^\ w^'^'^\ S,j) and four equations, three 
of which are linearly independent, namely three Einstein field equations (5-7) and the Bianchi 
identity. The system is thus initially undetermined and we need additional constraints to close 
the system. We have assumed that the DE is minimally interacting, T^'^''^ .^ — 0; thus due to 
the Bianchi identity (8) the perfect fluid component is also minimally interacting, T^™' .^ — 0. 
Hence, the Bianchi identity has been split into two separately additive conserved components; 
namely, the conservation of the energy-momentum tensor of the DE 

T(<i<=)"% = p(dc) + (l + ^(de)) ^(dc) [a ^ ^b\ ^ ^(,,) fA ^ 2^|\ ^ ^ ^^^^ 

and the conservation of the energy-momentum tensor of the perfect fluid component 



T^ 



V = pM + (l + ^(-))p(-)^ + 2| =0. (11) 



Once more we may split the conservation of the energy-momentum tensor of the dark-enery into 
two parts: 

y(dc)A"^^^ ^ r'(d^)^^, + r('i°)^% - 0, (12) 

where r'^^''^ -^ is the last term of the T^'^^^ ^^ in (10) and arises due to the deviations from 
'W^'^'^\ and T '^^''^ -i, is the the deviation-free part of the T^^''^ -^ in (10). Now, we will do the 
following strong assumption, 

r('^°^"V=p('^°^('4 + 27f)=0, (13) 

which also results in the deviation-free part of the T'-'^^^ .,j to be null, that is, 

y'(dc)A'-^^ = p(^^) + (l + T^f-^^)) pf'i^) (^ + 2^ = 0, (14) 



which looks hke the conservation of the energy momentum tensor of a minimally interacting 
perfect fluid. According to (13) and (14) the behaviour oi p^'^'^^ is controlled by the deviation-free 
part of the EoS parameter of the DE (w*^ '*")), but the deviations will affect p^'*") indirectly, since, 
as can be seen later, they affect the value of 'w^'^'^\ If the deviation parameters are assumed to 
be constants, to assure r''^'^-' -u = either S and 7 are trivially null or the ratio of the expansion 
rate on x axis to the expansion rate on y axis is equal to —2j/S which is a very special case. On 
the other hand, we may get more general solutions, if 5 and 7 are allowed to be function of the 
cosmic time t and we constrained S and 7 by assuming a special dynamic which is consistent 
with (13). The dynamic of the deviation parameter on the x axis, 5{t), is assumed to be 

^ ' 3B yA Bj p(d<=)' ^ '' 

and thus from (13) the deviation parameter on the y and z axes, 7(i), is found as 

lAfAB\l 
7(t) = -n-- - + 2— — -^. 16 

'^ -^ 3AyA B j p(d<=) ^ ' 

In such an assumption 5{t) and 7(t) are dimensionless parameters, and n is a dimensionless 
constant that parametrizes the amplitude of the deviation from w^^^' and can be given real 
values. The measure of the anisotropy of the DE may be given by {8{t) — ^{t))/w'''^°^ and it is 
null, which implies the DE is isotropic, when n — Q. 

The EoS parameter of the perfect fluid has been assumed to be constant, 

p(m) ^ ^ 

while w^'*®) is allowed to be a function of the cosmic time, since the current cosmological data 
from SNIa, CMB and large scale structures mildly favor dynamically evolving DE crossing the 
PDL as mentioned in the Section 1. Hence, since w'^'^'^\t) hasn't been constrained, we still need 
one more constraint to close the system. We imposed a law of variation for the Hubble parameter. 
The law of variation for the Hubble parameter which was initially proposed by Herman for RW 
space-time and yields a constant value of deceleration parameter [SH [5S1 [SHI [S71 [551 ES] • Such 
a law of variation for Hubble parameter given by Herman is not inconsistent with observations 
[551 I56j and is also approximately valid for slowly time varying deceleration parameter [57] . 
Recently Singh and Kumar proposed a similar law of variation for the Hubble parameter in 
anisotropic space-time metrics that yields a constant value of the deceleration parameter, and 
generated solutions for Bianchi Type-I [SHj , LRS Bianchi Type-II [SH [3S] , Bianchi Type-V [57] 
metrics in General Relativity. According to the proposed law, the variation of the mean Hubble 
parameter for the LRS Bianchi-I metric may be given by 

H = k{AB^)-'^, (18) 

where A: > and m > are constants. The spatial volume is given by 

V = a^ = AB^ (19) 



where a is the mean scale factor. The mean Hubble parameter H for LRS Bianchi-I metric may 
be given by 

i^-^--^--f4 + 2^V (20) 

a 3V 3yA Bj ^ ^ 

The directional Hubble parameters in the directions of x, y and z respectively may be defined 
as 

A B 

H^ = — and Hy=H, = —. (21) 

A y "■ B ^ ' 

(Since Hy — H^, in the following they are represented by Hy^z-) The volumetric deceleration 
parameter is 

9-- — - 22 

On integration, after equating (18) and (20), we get 

AB^ = cie^'^* for m = (23) 

and 

AB"^ = (mkt + C2) " for m^Q. (24) 

Here ci and C2 are positive constants of integration. Using (18) with (23) for to = 0, and with 
(24) for TO, ^ mean Hubble paremeters are 

H = k for m = (25) 

and 

H ^ k{mkt + 02)'^ for to y^ 0. (26) 

Using (23-24) and (19) in (22) we get constant values for the deceleration parameter for the 
mean scale factor as: 

g = TO, — 1 for TO, 7^ (27) 

and 

q = -1 for m == 0. (28) 

The sign of q indicates whether the model accelerate or not. The positive sign of q (i.e. to, > 1) 
corresponds to decelerating models whereas the negative sign — l<g<0for0<m<l 
indicates acceleration and q = for to, = 1 correponds to expansion with constant velocity. 

Using the deviation parameters (15) and (16), and the mean Hubble parameter (20) in the 
subtraction (6) from (7) we may get 

Tf4--Vf4--l3H^3ni7^ (29) 

dt yA B j yA Bj ^ ' 

after little manipulation. On integration of (29) by considering (25) and (26) we obtain 

4 - - = Ae-^''* + nk for to = 0, (30) 

A B ^ ' 
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A B X 3nk 

for m / and 3 (31) 



^ B (mkt + C2)^ (3 - m){mkt + C2) 

and 

A B A nk\n{3kt + C2) 



for m = 3, (32) 



^ B 3kt + C2 3fcf + C2 

where A is the real constant of integration. One can observe that A and n are two parameters 
that parametrize the difference between the directional Hubble parameters. Now, we can find 
A{t) and B{t) explicitly for all m values by using (30-32). 

3.1 Model for m = {q = -1) 

Using (30) we may get the ratios of the scale factors A{t)/B{t), and on manipulating the result 
by using (23) we get the following exact expressions for the scale factors: 

Bit) = cii/3K-i/3efc*+H«-'^'''-§«fe*, (34) 

where k is the positive constant of integration. The spatial volume of the universe is found as 

V = cie^*^*. (35) 

The directional Hubble parameters as defined in (21) are found as 

H^ = k+ -Ae-3fc* + -kn, (36) 

3 3 

Hy.^k- -Xe-^''* - -kn. (37) 

y,2 3 3 V ; 

The anisotropy parameter of the expansion A is defined as 

3 / TT TT\ 2 



^iEHt^' (38) 



A- , , 

3^ V H 

where Hi (i=l,2,3) represents the directional Hubble parameters in the directions of x, y and z 
respectively. By using equations (25), (36) and (37) in (38) we get 

A =9^ ^^ (39) 

for the anisotropy of the expansion. The expansion scalar, defined by = u^-j, is found as 

e = 3fc = 3H. (40) 

The shear scalar, defined by ct^ = ho'ijcr^-' where aij ~ Ui-j + Uj-^i — ^giju'^-k is the shear tensor, 
is found as 

a2 = i(Ae-3fe* + nfc)'. (41) 



Using the scale factors in (11), the energy density of the perfect fluid is found as 

pM(t)=.p('")e-3'=(i+-'"')*. (42) 

The energy density of the DE may be found from (5) by using the scale factors and the energy 
density of the perfect fluid (42), and may be written in terms of A{t) as 

p('i°)(t) = 3fc2(^l-iA(t))-pM(t). (43) 

Using the scale factors and (43) in (14) we get 

, , ^ , , 3w('^V^"'Hi) + A^e-efc* + fc2 (9 - n^) , , 

^' 3p('°)(t) + A2e-6fc* + 2Afcne-3fct_A:2(9-n2) ^ ^ 

for the deviation-free EoS parameter of the DE. Using the scale factors and (43) in (15) and 
(16) we may get deviation parameters as following, 

2Afcne-3fc* + 2k^n[n - 3) 

^ ' ~ 3p('°)(t) + A2e-6fe* + 2Afcne-3'=^* - fc2 (9 - n2) ' ^ ' 

2Afcne-3'=* + 2fc2n(n + 3/2) 
^^' ~ 3/9('")(i) + A2e-6'=* + 2Afcne-3'=* ~ k'^ {% - n"^) ' ^ ' 

3.2 Physical behaviour of the model for m = (g = — 1) 

For this model q = —1 and dH/dt = 0, which implies the greatest value of the Hubble parameter 
and the fastest rate expansion of the universe. Thus, this model may represent the inflationary 
era in the early universe and the very late times of the universe. 

The spatial volume V is finite at i = 0, expands exponentially as t increases and becomes 
infinitely large at i = cxd. The directional Hubble parameters H^ and Hy^z are finite at i = 
and t ~ oo. They deviate from the mean Hubble parameter H due to A and n. While A is 
supporting (opposing) the expansion on the x axis, it opposes (supports) the expansion on y 
and z axes. However it loses its effect exponentially by the cosmic time t. The anisotropy of 
the DE has a similar effect on the directional Hubble parameters, but its effect persists. While 
n is supporting (opposing) the expansion on the x axis, it opposes (supports) the expansion on 
y and z axes. The anisotropy of the DE doesn't always act so as to increase the anisotropy of 
the expansion. In fact, when the signs of A and n are opposite the overall effect is to lower the 
expansion anisotropy. 

The expansion scalar is constant throughout the evolution of the universe. The shear scalar 
is also finite at i = and tends to n^k'^/3 as t increases. If —3/2 < n < 3 all the axes will 
expand to infinitly large values as t — ^ cxd. On the other hand, the space-time exhibits a pancake 
type singularity for n < —3/2 and a cigar type singularity for n > 3 at t = oo. If n = —3/2 
while the x axis converges to a constant, y and z axes expand to infinitly large values, and it is 
vice versa if n = 3. 

The energy density of the perfect fluid p^™' decreases exponentialy and converges to zero 
since w^™' > by definition, while p*-'*°^ changes slightly at early times and converges to a non- 
zero value as t increases. Thus, the ratio of /9'^'^°)/(p('"^ -fp'-'*"-') converges to 1 as i increases, that 
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Figure 1: p'^'^°^ versus cosmic time t for different Figure 2: w^'^") versus cosmic time t for differ- 
values of n in model ?ti = 0. A and k have been ent values of n in model tti = 0. A and k have 
chosen as 1. been chosen as 1. 



is the DE dominates the perfect fluid in the inflationary era as expected. Since DE dominates the 
perfect fluid in the inflationary era as mentioned above, we may neglect the perfect fluid while 
examining the properties of the DE in this model. p^'^°'^ is always positive provided |n + A/fc| < 3 
and \n\ < 3, begins with a finite value at i = 0, exhibits different behaviours depending on the 
choice of the parameters at the earlier times and tends to 2>k^{l — ri^/9) for large values of t; see 
Fig[TJ It can be observed that the bigger the anisotropy of the expansion the lower the p'--'^^^ in 
any given instant. A loses its effect exponentially as t increases, thus we may say that the higher 
the \n\ the lower the p^^^' in any given instant for relatively big t values. The EoS parameter of 
the DE w'^°) may begin in phantom {w < —1) or quintessence (w > —1) region and tends to —1 
(cosmological constant, w = —1) by exhibiting various patterns as t increases; see Fig [2] One 



can observe that p^'^'^'' increases when 



n(do) 



< 



-1, decreases when w^'^° 



> 



-1 and is constant 



when w'-^^^ = —1, as would be expected. 

The anisotropy of the expansion decreases monotonically as t increases when n — 0. However, 
it exhibits nontrivial behaviour at the early times of the universe and converges to a non-zero 
value for the late times when n ^ 0; see Fig[3l The deviation parameters S and 7 are finite at 
t = 0, and converge to 2n/(n-f 3) and 7i(2ri-|-3)/(7i^ — 9) respectively as i ^ 00. The anisotropy 
of the DE which has been defined as {S — 7)/w'^'^°^ tends to 9n/(n^ — 9) as t increases; see FigH) 



3.3 Model for m ^ {q ^ -1) 

The solutions in this subsection are valid for all possible values of m except for m = 3 and 
m — 0, thus the solutions for m ~ 3 arc given in the following subsection. 

From (24) one can see that the initial time of the universe is i* = —C2/mk for m ^ 0. For 
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brevity of the equations, we may redefine the cosmic time as 

t' = mkt + C2, 



(47) 



and by doing that the initial time of the universe has also been set to t' = 0. Thus we may 
rewrite the metric as 



ds^ = {mk)-^dt'^ - A{t'fdx^ - B{t'f{dy^ + dz^) 



(48) 



Using (31) we may obtain the ratios of the scale factors A{t)/B{t), and on manipulating the 
result by using (24) we get the following exact expressions for the scale fators: 



T/o ,1 2re 2 A f'l-- 



A(i') = K^i''" '"("-3)e3M^^^^^ 



(49) 



T /o ,_LJ li 1 A ../l-T^ 



B(t') = K^'i' ""<"•"'> e" 3 M^^r^^ "^ (50) 

where k is the positive constant of integration. The spatial volume of the universe is found as 



The directional Hubble parameters as defined in (21) are found as 

w !,+'-! _u 2 \+'-^ 2nfc 1 



(51) 



(52) 
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Hy^. = kt'-' \\t'-^ + j^,t'-\ (53) 

Using (26), (52) and (53) in (38) we get 

for the anisotropy parameter of the expansion. The expansion and shear scalars are, respectively, 
found as 

e = 'ikt'^^ = 3H, (55) 

...>A-(^,.-*__!^)^ „„ 

Using the scale factors (49-50) in (11), the energy density of the perfect fluid is found as 

p(-)(i')=p('")i'-™<'+"'*"''. (57) 

The energy density of the DE can be found from (5) by using the scale factors (49-50) and the 
energy density of the perfet fluid (57), and may be written in terms of A(t) as 



p(^^)(t') = 3fc2 fl - -Ait')\ t'"' - p^"'Ht'). 
Using (49-50) and (58) in (14) we get 



(58) 



w^'^Ht') = ^^^ 7 T T^-^ \ ^^^^^ -— -, (59) 

for the deviation-free EoS parameter of the DE. And finally using equations (49-50) and (58) in 
(15) and (16) we may get deviation parameters as following, 

-2 fou 2 \+/l--^ , 2nk 



nkt'^"- 2k- iXt' 



3 ' (m-3) 

3fc2 (l - iA(t')) t'"^ - M'"H^') 



m - .,..\. :,.„-. ::z:^ (60) 



n. 



kt'-^-k-^Xt''-^ ' 2„fc 



/,/N ^ \ 3 (™-3) / , , 

^^^" 3fc2 (1 - iA(t')) t'"' - ^('"Hi') ■ 

3.4 Physical behaviour of the model for m ^ {q ^ —1) 

The universe accelerates for < to < 1, decelerates for m > 1 and expands with constant 
velocity for to, = 1. This model may represent the radiation dominated era for to = 2 and 
y;(m) _ i^^^ gj^^ |^]-^g matter dominated era for to = 3/2 and w^™) — 0. 

The mean Hubble parameter H, expansion scalar 9 and shear scalar a^ are infinitly large at 
t' = 0, and null at t' = oo. 
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If m > 3, the space time exhibits a cigar type singularity at i' = cx) for A > and a pancake 
type singularity for A < 0. If m < 3, all the axes expand to infinitly large values provided 
(m — 3)/2 < n < 3 — m as t' ^ oo. When n = 3 — m, A{t') takes infinitly large values as 
t' -^ cxo, while B{t') converges to k~^/'^. When n — (m — 3)/2, B{t') takes infinitly large values 
as t' -^ oo, while A{t') converges to k^^^. When n > 3 — m space time exhibits a cigar type 
singularity at t' — oo. When n < (m — 3)/2 space time exhibits a pancake type singularity at 
t' = oo. 

A diverges as t' -^ 0, converges to a constant as i' — > oo for m < 3 and vice versa for 
m > 3. One can observe that the anisotropy of the expansion lowers the value of p^'^®^ . Thus, A 
must be considered while examining the dynamics of the energy density of the DE, particularly 
since A diverges in the above mentioned limits. In the model for m < 3, p'-'^''^ will eventually 
obtain negative values as t' goes to zero due to the divergence of A. Hence, this model is not 
appropriate for representing the relatively earlier times of the universe. On the other hand, this 
model can set to represent relatively later times of the universe since we can always set p^'^'^'> 
to be p^'^"'' > as t' -^ oo by choosing the parameters in suitable intervals. Similarly, in the 
model for ttt. > 3, p^'^") will eventually obtain negative values as t' -^ oo. Hence, this model is 
not appropriate for representing the late times of the universe, but for the earlier times of the 
universe by choosing the parameters in suitable intervals. 

w^'^®) and {S — "f)/w'^'^'^\ say the anisotropy of the DE, exhibits various dynamics according 
to the choice of the parameters. It is worth to mention here once more that the anisotropy of 
the DE doesn't always act so as to increase the anisotropy of the expansion, when the signs of 
A and n are opposite the overall effect is to lower the expansion anisotropy. 

We may try to set t' = 1 to the present universe by using the recent cosmological data which 
favor Hq « 0.71 (dimensionless Hubble parameter), fig °^ « 0.73 (density parameter of the dark 
energy), fig = 1 — fig « 0.27, Wq ss —1 ^5D|, w^"^^ — and a space-time that can be 
represented by flat RW metric which implies Aq — for the present universe, where o subscripts 
represent present (t' = 1) values of the parameters. With above parameters we may obtain 

..(c) 

present critical energy density of the universe as pg — 1.5123, the deceleration parameter 
as go = -0.595 which implies m = 0.405 in our model, p^,'"^ = 0.408321, p^,'^"^ = 1.103979, 
^{dc) = _i and /^jj ^ by choosing A = -0.8208092487?i. 

In this model with the above values of the parameters, which may correspond to the present 
universe in the vicinity of i' ~ 1, p^™> is roughly proportional with t'^^, while p^'^''' is roughly 
proportional with t~^ and gets its highest value at i' = 1 where A = 0. That is, the DE 
eventually dominates the perfect fluid. Some graphs may be plotted with above values, since 
we can think that q is changing only slightly for the present universe, for n = 0, n — —0.3 
and n = 0.5. Plots has been extended towards to the later times more than towards to the 
earlier times since, as mentioned above, models with m < 3 are not appropriate for representing 
relatively earlier times of the universe, p^'^'^' increases when w^^^' is in the phantom region, 
begins to decrease when w^^^' passes into the quintessence region and converges to zero as 
t' -^ oo; see Fig [5] w^^°^ begins in the phantom region, increases and becomes —1 at i' = 1 
then passes into the quintessence region and tends to a constant which is in the quintessence 
region; see Fig [51 p^™-* decreases as t' increases independent of n; see Fig[71 Q,'^'^^^ increases as 
t' increases, becomes 0.73 at i' = 1 and converges to 1 as t' keeps on increasing; see FiglH A 
decreases as t' increases and becomes null at i' = 1 then tends to a constant which depends 
on n, see Fig |9l The bigger the A the lower the p*-'^°\ consequently lower the fi^'^") in any 



13 





-n = 



n = -0.3 



- n = 0.5 



-n = 



n = -0.3 



- n = 0.5 



Figure 6: The plot of u'^'^''-' versus cosmic time 
Figure 5: p'^'^°^ versus cosmic time t' in the t' in the vicinity of i' = 1 in the model m = 
vicinity of i' = 1 in the model m — 0.405 for 0.405 for different values of n. A has been set to 
different values of n. A has been set to nuU at null att' = 1 by choosing A ~ -0.8208092487n. 
t' = 1 by choosing A = -0.8208092487n. w'-'^^'^l) = -1 

given instant. The anisotropy of the DE, {S — ^)/'w'^'^°\ is null for n — 0, changes slightly and 
converges to a non-zero constant for n ^ 0; see FigfTOl 



3.5 Model for m = 3 (q = 2) 

From (24) we can see that the initial time of the universe is t, ~ — C2/3fc for rn — 3. For brevity 
of the equations, we may redefine the cosmic time as 



t' ^3kt + C2, 



(62) 



and by doing that the initial time of the universe has also been set as t' = 0. Thus we may 
rewrite the metric as 



ds^ = {3k)-^dt'^ - A{t'fdx^ - B{t'f{dy^ + dz^). 



(63) 



Using (32) we may get the ratios of the scale factors A{t)/B{t) for m = 3, and manipulating 
the result by using (24) we get the following exact expressions for the scale fators; 



B{t') =:K-h'*^^*e~T5'"(*')', 



(64) 
(65) 
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Figure 8: The plot of Q^^'^'' versus cosmic time 
Figure 7: The plot of p^™^ with w^™) = 0, vcr- t' in the vicinity of i' = 1 in the model m — 
sus cosmic time t' in the vicinity of t' = 1 in 0.405 for different values of n. A has been set to 
the model m = 0.405. The behaviour of the null at f = 1 by choosing A = -0.8208092487n. 
p^"") is indepedent of n. rj('i°)(l) = 0.73 

where k is the positive constant of integration. The spatial volume of the universe is found as 

V = t'. (66) 

The directional Hubble parameters as defined in (21) are found as 



H^^kt' ' + -{X + nk\n{t'))t' \ 
Hy^, = kt'^^ - ^ (A + nk\n{t')) t''\ 



Using (26) and the scale factors in (38) we get 



A= -k'^ {nk\n{t') + Xf 
9 



(67) 
(68) 

(69) 



for the anisotropy parameter of expansion. The expansion and shear scalars are, respectively, 
found as 

9 = 3kt'^^ = 3H, (70) 

^ '- '^'■'"' (71) 



a^ = -(nfcln(t') + A) t'" . 
Using the scale factors in (11), the energy density of the perfect fluid is found as 



(72) 
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Figure 10: Plot of the anisotropy of the DE 
Figure 9: The plot of A versus cosmic time t' versus cosmic time t' in the vicinity of t' = 1 
in the vicinity of t' = 1 in the model m = 0.405 in the model m = 0.405 for different values of 
for different values of n. A has been set to null n. A has been set to null at t' = 1 by choosing 
at i' = 1 by choosing A = -0.8208092487n. A = -0.8208092487n. 



The energy density of the DE can be found from (5) by using the scale factors and the energy 
density of the perfet fluid (72), 



p('i°)(t') = 3fc2 (^1 - iA(t')) t'^' - p('"Ht'), 



(73) 



Using (64-65) and (73) in (14) we get 

(do). ,N _ 3fc2t'-2 _ 1 (nfcln(t') + A)2f'-2 + lnk(nk\xi (f) + A)t'"2 - u;('»)p('")(t) 

"" ^* ^ " 3fc2 (1 - iA(i')) t'"' - p('"Hi') 

for the deviation-free EoS parameter of the DE. And finally using equations (64-65) and (73) in 
(15) and (16) we may get deviation parameters as following, 



(74) 



m 



lit') = 



nfc(2fc-|nfcln(i')-|A)i'"^ 
3fc2(l-iA(i'))t'"^-M'"H^')' 
nk {-k - InkHt') - l\) t'-^ 
3fc2(l-iA(t'))t'"^-/o('"Hi')' 



(75) 
(76) 



3.6 Physical behaviour of the model for m = 3 (g = 2) 

Nothing physically special has been observed in this model for further investigation. But we 
may mention that this model may only be valid for intermediate epochs of the universe. Because 
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A diverges at i' = and i' = oo in this model, thus p^'^®^ will eventually get negative values as 
t' goes to extrem values, since the anisotropy of the expansion contribute the energy density of 
the DE negatively. 

4 Conclusion 

Locally rotationally symmetric Bianchi-I cosmological models with dynamically anisotropic dark 
energy (DE) and perfect fluid have been constructed in General Relativity. We assume that 
the dark energy (DE) is minimally interacting, has dynamical energy density and anisotropic 
equation of state parameter (EoS). The conservation of the energy- momentum tensor of the 
DE has been assumed to consist of two separately additive conserved parts. A special law has 
been assumed for the deviation from isotropic EoS, which is consistent with the assumption 
on the conservation of the energy-momentum tensor of the DE. Exact solutions of Einstein's 
field equations have been obtained by assuming a special law of variation for the mean Hubble 
parameter, which yields a constant value of the deceleration parameter and is not inconsistent 
with observations. Some basic geometrical and kinematical features of the models and the 
dynamics of the anisotropic DE in these models have been examined. 

A more general EoS parameter has been introduced for the DE, and the isotropic DE can 
be recovered by choosing n to be null, where n parametrizes the amplitude of the deviation 
from isotropic EoS parameter of the DE. In all the models, while the anisotropy of the DE 
contributes to the expansion of one of the scale factors, it opposes to the expansion of the other. 
Two parameters that parametrize the difference between directional Hubble parameters, A and 
n, give rise to non-trivial, dynamically anisotropic, expansion histories. This allows for the 
possibility to fine tune the isotropization of the Bianchi metric during an accelerating epoch of 
the universe in order to generate arbitrary ellipsoidality by choosing a suitable value of n. Such 
a result provides also the possibility to fine tune the CMB anisotropy [571 HH] ■ The anisotropy 
of the expansion can mildly or totally isotropize in relatively earlier times of the universe. 
Nevertheless, it converges to a nonzero constant value for the later times of the universe in all 
models. It is also observed that the anisotropy of the DE energy doesn't always act so as to 
increase the anisotropy of the expansion, when signs of A and n are opposite the overall effect 
is to lower the expansion anisotropy in relatively earlier times of the universe. 

One interesting observation is that the higher the anisotropy of the expansion the lower the 
energy density of the DE in each given instant of the cosmic time. u'''^°-' is also dynamical and 
exhibits non-trivial behaviour in our model, but it eventually converges to a constant in the 
quintessence region for m ^ and converges to —1 for m — Q. 

Such cosmological models are of interest because they give rise to an ellipsoidality of the 
universe in spite of the inflation [321 SSI HD], which is one of the promising proposals to the 
solution of the quadrupole problem [5T1 [521 HZl HH] and can also be checked by the direction 
dependency of the redshift-luminosity relation of the SNIa observations [49l |50| . 
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